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ABSTRACT 

Recent numerical simulations of binary black hole mergers show the possibility of producing very 
large recoil velocities (> 3000 km s _1 ). Kicks of this magnitude should be sufficient to eject the final 
black hole from virtually any galactic potential. This result has been seen as a potential contradiction 
with observations of supermassive black holes residing in the centers of most galaxies in the local 
universe. Using an extremely simplified merger tree model, we show that, even in the limit of very large 
ejection probability, after a small number of merger generations there should still be an appreciable 
fraction (> 50%) of galaxies with supermassive black holes today. We go on to argue that the inclusion 
of more realistic physics ingredients in the merger model should systematically increase this retention 
fraction, helping to resolve a potential conflict between theory and observation. Lastly, we develop 
a more realistic Monte Carlo model to confirm the qualitative arguments and estimate occupation 
fractions as a function of the central galactic velocity dispersion. 
Subject headings: black hole physics - galaxies: nuclei 



Or 

6 



> 
oo 

o 
i> 
o 



X 



1. INTRODUCTION 

In the standard model of hierarchical structure forma- 
tion, the most over-dense regions of the universe collapse 
at early times to form small gravitationally bound clumps 
of dark matter and baryons. These proto-galaxies sub- 
sequently merge with each other and form larger and 
larger objects, up to the galaxies and clusters we see 
today dPress fe SchechTer1ll974t lEfstathiou fc Reeslll988t 
lHaehnelt fc Reesl 119931 : ICole et alJl2000j ). Observations 
of the local universe suggest that most galaxies at red- 
shift z = are also host to a supermassive black hole 
(SMBHs) near their centers (jRichstone et al.lll998f ). Ob- 
servations of distant quasars in the Sloan survey sug- 
gest that th ese SMBHs are common up to redshifts of at 
least z ~ 6 ()Fan et al.l l2001). One natural explanation is 
that the SMBHs originally formed in the early universe 
and then simply followed their host galaxies through- 
out multiple generations of hierarchical mergers: when- 
ever a pair of galaxies merged, their central black holes 
would sink together through dynamical friction and even- 
tually merge themselves, resulting in a single, more mas- 
sive galaxy with a single, mor e massive BH at its center 
(jKauffmann fc Haehneltll2000HCattaneoll2001[ ). 

In the absence of special symmetries, the final black 
hole will receive a linear momentum recoil from the gravi- 
tation al radiation emit t ed during the b inary merger pro- 
cess (jBekensteinl 1 1 9 73t iFitchett 1 1 1 9831 ) . Recent numer- 
ical relativity simulations suggest that this recoil can- 
not exceed 175 km s ~ 1 for no n- spinn i ng black holes 
(iHerrmann et al.ll2006HBaker et al.ll2006tlGonzalez et all 
l20n(t) . but may be as large as ~ 3600 km s 1 for 
maximally spinning equal-mass b lack holes with spins 
anti-aligned in the orb ital plane (|Gonzalez et alj l2007t 
iCampanelli et all [2007). In this latter case, the re- 
sulting BH would almost certainly be ejected from the 
gravitational well o f its host galaxy (|Menou etall 12001 : 
iMerritt et aTll2004ft . 

However, since observations of the local universe sug- 
gest that SMBHs are quite common, there appears to 
be something preventing these ejections. As a first at- 



tempt at resolving this conflict. ISchnittman fc Buonannol 
(|2007f ) carried out Monte Carlo simulations of the bi- 
nary BH recoil for a large range of mass ratios and spin 
orientations, finding a relatively small fraction /iooo — 
0.05 — 0.1 of systems with kick velocities > 1000 km s _1 . 
These fractions correspond to black holes with dimen- 
sionless spin parameters of a/M = 0. 9, which appear 
to be typical from AGN observations (| Yu fc Tremaind 
| 2002HElvis et alJl2002tlWang et al.ll2006h . More recently, 
iBogdanovic et all (|2007f ) proposed a mechanism in which 
a circumbinary gas accretion disk would align the spins 
of the two black holes with the orbital axis of the disk 
prior to merge r, thereby limiting the maximum recoil to 
< 200 km s" 1 (jBaker et al.ll2007h . 

In this paper, we take a rather different approach to 
the problem and question the very premise of the appar- 
ent conflict: Does a large recoil velocity necessarily imply 
a small occupation fraction? Furthermore, does a small 
occupation fraction at high redshift imply that only a 
small number of SMBHs will survive until today? If, for 
example, 95% of all BH mergers result in an ejection from 
the host galaxy, does this really mean that at least 95% 
of galaxies today should be devoid of SMBHs? Using a 
binary tree merger model for the hierarchical growth of 
BHs and their host galaxies, we show that the answer to 
these questions is a resounding "NO!" Our simple ana- 
lyt ic predict i ons a ppear to agree well with recent work 
bv lVolonteril (|2007f l. which employs a much more detailed 
merger evolution, including BH evolution due to accre- 
tion, but also finds a large occupation fraction at low 
redshift. 

The formal mathematical argument is given below in 
Section [2J but a simple qualitative reasoning is this: in 
every generation of the merger tree, the number of galax- 
ies decreases, so the fraction with BHs can easily increase. 
If every BH host galaxy merges with an empty galaxy, 
then the fraction of galaxies with SMBHs can double 
from 0.5 to 1. Including the losses due to ejected BHs, 
we find a steady state solution with occupation fraction 
/ = 1/(1 +p e j), where p Ci is the probability that a merg- 
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ing BH is ejected due to gravitational recoil. In Section 
we discuss a number of potential physics ingredients 
that could be added to the simple model, arguing that 
each one would systematically increase the fraction of ob- 
served SMBHs today. In Section[4]we attempt to confirm 
some of these arguments with a more astrophysically re- 
alistic Monte Carlo merger model, and in Section [5] we 
present our conclusions. 



2. SIMPLE MERGER TREE MODEL 

We begin with the simplest possible model for a galac- 
tic merger tree: a binary tree where in each generation 
every galaxy merges with one "spouse" and produces a 
single "child" (see Figs. Q] and [2] for a schematic). Each 
parent galaxy may or may not have a SMBH at its center, 
and if they both do, then the respective black holes are 
assumed to merge into a single black hole, which is then 
ejected due to gravitational recoil with a probability of 
p c j. In the zeroth generation, before any mergers, every 
galaxy is assumed to have a central black hole present. 
Thus in the first generation, a fraction fx = 1 — p C j of 
galaxies will have SMBHs. 

For a galaxy in the (i + l) st generation, the probability 
of having a black hole is 

/,+i=0x(l-/,) 2 

+/i(l-/i) + (l-/i)/< 

+ (l-Pej)/f 
= /i(2-/ i - Pej / i ), (1) 

where the first line corresponds to having two parents 
without black holes, the second line represents either a 
black hole mother or father, and the third line is the 
probability of both parents having black holes, but the 
child black hole not getting ejected in the merger. Taking 
the convergence limit of fi+\ = fi for large i, we find 



This remarkable result suggests that, in the limit of 
infinite recoil velocity, 50% of all galaxies will still retain 
a central black hole. This is possible because, during 
one generation of galaxy mergers with /j = 0.5, while 
effectively half of the black holes are removed by the 
recoil, the total number of galaxies is also cut in half, 
thus maintaining a constant fraction of SMBHs. 

In Figures [1] and [2] we show examples of the binary 
merger trees that obey the simple rules described above. 
The heavy black circles represent black holes, while the 
empty circles are galaxies without SMBHs. For ejection 
fractions p c) = 0.1 (Fig. CD) and p c ; = 0.9 (Fig. [2|) we 
show six generations of galaxy mergers, beginning with 
/i = 1 — p cy In both cases, the system approaches the 
equilibrium fraction f^ of equation ([2|) after only a few 
generations. 

For a more quantitative estimate of the "convergence 
time" that it takes to reach equilibrium, we see from 



equation (CQ) that 
/i = l -p C y 
h = 1 - Pcj + P% ~P%, 

/ 3 = 1 -pej + P% -Pcj +Ptj -P% +P% -Pip 

f. £< 1 ?'/<:• (3) 

3=0 

This series not only converges to equation ([2]) as ex- 
pected, but does so at quite a rapid pace. Even for ex- 
tremely large ejection rates of p c j — 0.95 (i.e. /i = 0.05), 
fi is within 5% of its asymptotic value of /oo = 0.513 
after just five generations. Because of this rapid conver- 
gence, we find that the final results are independent of 
the initial occupation fraction /q. 

3. ADDITIONAL PHYSICS COMPONENTS 

The model as described above is based on two major 
simplifying assumptions: a perfect binary merger tree, 
and a constant ejection probability. Of course, in reality 
both of these assumptions are in all likelihood invalid. 
Here we list a number of possible modifications based on 
more realistic astrophysics, and for each one argue that 
their inclusion will only increase the fraction of SMBHs 
observed today. 

The binary merger tree assumes that for a given gen- 
eration, every galaxy has exactly the same mass and 
merges with exactly one other galaxy, halving the to- 
tal number of galaxies in each subsequent generation. In 
practice, most hierarchical merger simulations go back- 
wards i n time, using the extend ed Press-Schcchter for- 
malism ()Press &: Schechterll 19741) to estimate the progen - 
itor masses at a given redshift (jVolonteri et alJ [2003). 
Many of these simulations find that a single trunk of 
the merger tree dominates, with a large n umber of small 
branc hes joining in at different redshifts (jMalbon et al.l 
2006). This means that the mass ratio for a typical 
(proto-)galaxy merger can be significantly different than 
unity. Assuming that the SMBH masses scale with 
their host mas ses, this suggests th at BH mass ratio will 
not be unity (jSesana et al.ll2007l ). While the Fitchett 
scaling for non-spinning black holes favors a mass ra- 
tio o f about three-to-one for maximum recoil (|Fitchett I 
1983), numerical simulations imply that the super- 
massive kicks of j> 1000 km s -1 requ i re nearly-equal 
masses and spins (iHerrmann et al.ll2007t iGonzalez et al.l 
l2007HCampanelli et alJl2007L iTichy k. Marronettill2007h . 
so a wider range of mass ratios will tend to decrease the 
ejection probability p G j . 

Furthermore, as the galaxies evolve along the hierar- 
chical merger tree, their masses will increase, and thus 
so will their escape velocities. Of course, the BHs are 
also growing in mass through mergers and accretion dur- 
ing this time, but the recoil problem is strictly scale- 
invariant: the kick velocity is a function only of the mass 
ratio and the dimcnsionlcss spin parameters. Thus if the 
masses of the BHs and their host bulges all double, the 
recoil will be the same, but the escape velocity will be 
larger, and thus p Ci will be smaller. In this scenario, even 
if most of the BHs are ejected at large redshifts, after a 
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few generations /, will grow rapidly due to a large num- 
ber of "single-parent" mergers, ultimately converging to 
a larger population fraction based on the smaller value 
of p j at late times. 

Another effect that may increase the escape velocity 
is the formation of a triple-BH system by the relatively 
prompt merger of three h aloes before their respe ctive 
BHs have time to merge (iHoffman fc Loebl I200TO . In 
many of these cases, one of BHs is ejected by Newto- 
nian three-body interactions, leaving a "normal" binary 
BH in a more massive host galaxy (i.e. made of three 
parents instead of two), further increasing the escape ve- 
locity and lowering the ejection probability. However, 
this same three-body interaction may be strong enough 
to give both the single and binary a large enough New- 
tonian recoil so that all three BHs are ejected. 

In addition to these modifications to the merger tree 
physics and escape velocities, there are also a number 
of processes that more directly affect the actual recoil 
velocity. Since the largest kicks are found in systems with 
the BH spins in the orbital plane, any systematic effect 
that tends to avoid this orientation will thus reduce the 
expected kick velocity. One particularly strong influence 
on the spin orientation is the torque produced by a single 
circumbinary accretion disk, which can align both black 
hole spins with the orbital angular momentum with high 
efficiency (Bogdan ovic et aDl2007t ). In this orientation, 
the maximum recoil sh ould not be more than ~ 200 km 
s _1 ([Baker et al.ll2007l ). well below the escape velocities 
of most present-day galaxies. In the event that there is 
no surrounding accretion disk (i.e. a "dry merger"), the 
two black h ole spins may beco me aligned via spin-orbit 
resonances (ISchnittman I [200^ ). but this process likely 
requires somewhat special initial conditions. 

Lastly, there is also the possibility of creating a new 
SMBH ex nihilo during the galactic mergers, which typi- 
cally are accompanied by m assive gas inflows to t he cen - 
ter of the resulting galaxy iMihos &: Hernquistl (|1994| ). 
This rapid increase in gas density will trigger a burst 
of massive star formation, which may then proceed 
to form super-mass ive stars through runaway mergers 
(|Gurkan et al.ll2006t ) , in turn collapsing to form the seeds 
of SMBHs, which will also be surrounded by copious 
amounts of fuel to accrete more mass. In this way, or- 
phan black holes can appear in the merger tree, further 
increasing the overall occupation fraction. 

4. ASTROPHYSICAL MERGER TREE 

In an attempt to verify some of these qualitative 
claims, we now develop a slightly more physical Monte 
Carlo merger model, motivated by basic astrophysical 
arguments and observations of local SMBHs. Instead of 
using a constant ejection probability and a perfectly bi- 
nary merger tree of equal-mass BHs and galaxies, we now 
consider an initially flat BH mass distribution function 
dN/dM, ~ M^ 1 , and assume the galaxy mass is pro- 
portional to BH mass. If an occupied galaxy merges 
with an empty galaxy, we adjust the final BH mass 
to "agree" with the total mass of both parent galax- 
ies (physically represented by some accretion episode). 
For an isothermal sphere model, it turns out the galaxy 
mass is not important in determining the escape ve- 
locity, which depends only on the velocity dispersion: 
v csc — 2<t v . The velocity dispersion is in turn deter- 



mined from the BH mas s via the M — a relation of 
iFerrarese k Merrittl (|2000f ): iMerritt fc Ferraresel (|2001f ): 

/ v 4.7 

M. w 1.3 x 10 s ( — T ) Mq . (4) 

V200 kms" 1 / 

In each generation, a fraction f m of all the galaxies 
(selected randomly) merge. In the limit of f m — > 1 (as in 
the binary merger tree), the low-mass tail of the distribu- 
tion is prematurely depleted, and if f m — ► 0, no evolution 
takes place at all. However, we find that for the interme- 
diate range 0.25 < f m < 0.75, the net results are largely 
independent of f m . For each merger, we determine the 
final occupation as in the binary model: two empty par- 
ents create an empty child, a single BH mother or father 
will create a BH child, and in the case of two BH parents, 
the child black hole will be ejected if the kick velocity is 
greater than the escape velocity of the child galaxy. 

To determine the kick velocity in this more astrophys- 
ically rea listic model, we employ the an alytic fits pre- 
sented in ISchnittman &; Buonannol (|2007t ). assuming all 
BHs are rapidly spinning with a/M = 0.9. In this case, 
for a given mass ratio q = mXjml < 1, 90% of the merg- 
ers should produce kick velocities less than vqq: 

v 90 (q) w 7^— |J V(l - q) 2 + 1-4(1 + q? km s"\ (5) 

with the actual recoil Ukick selected randomly 
from the cumulative distribution function 
(|Schnittman fc Buonannoll2007t ) 

Pcdf(Wkick;g) = 10-^/^^. (6) 

In the upper panel of Figure O we plot the occupa- 
tion fraction as a function of a v , employing the model 
described above for two cases: an initial occupation of 
fo — 0.99 (solid curves) and for / = 0.4 (dashed curves). 
In both plots, the fractions /j refer to the average num- 
ber of mergers for each galaxy, not the number of gener- 
ations (Aggers » A gcn / m /(2 - f m )). After a few merg- 
ers, both cases converge to the steady-state case of /oo 
(heavy black curve), just as in the simple binary merger 
model. In the lower panel, we repeat the same calcu- 
lation, now multiplying wgo by a factor of three, to see 
the effect of extremely large kicks. While the occupation 
fractions clearly decrease somewhat, they remain signifi- 
cantly above 50%, just as predicted by the binary merger 
model. 

We do not claim that this improved model should 
be seen as reliable for making quantitative astrophysi- 
cal predictions, but rather present it as an application 
of the simple binary tree model and a confirmation of 
many of the qualitative arguments from Section O At 
the same time, it is noteworthy that the occupation frac- 
tions for the standard model (upper panel of Fig. [3]) seem 
to agree quite close ly with the much more detailed model 
of I Volonteril (|2007t ). We anticipate that future observa- 
tions should be able to measure this distribution at larger 
redshifts and smaller o~ v , ultimately explaining how the 
BH occupation fraction evolves in time. 

5. CONCLUSIONS 

Based on a very simple binary merger tree model, we 
find that the fraction of galaxies hosting a SMBH should 
be > 50% even in the limiting case of very large recoil 
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velocities (or alternatively, very small escape velocities). 
This would be particularly important at high redshifts, 
when typical galaxies have small escape velocities and the 
seeds of today's SMBHs are presumably formed. Includ- 
ing qualitative arguments about the hierarchical growth 
of the host galaxies, the ejection probability p e j will tend 
to decrease with cosmological redshift, further enhancing 
the fraction of central BHs observed today. Including 
spin-alignment effects from accretion will increase this 
fraction even more. 

In Section [4j we presented a somewhat more physical 
Monte Carlo model and were able both to confirm the 



predictions of the binary merger tree and also reproduce 
the ba sic results of the detailed calculations of IVolonterl 
(2007). We thus conclude that the very large kicks pre- 
dicted by numerical relativity should not necessarily be 
seen as contradictory to the seemingly ubiquitous popu- 
lation of galaxies with SMBHs observed today. 

The author would like to thank Alessandra Buonanno, 
Julian Krolik, and Cole Miller for helpful discussions and 
comments. Partial support comes from NSF grant PHY- 
0603762. 
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Fig. 2. — Merger tree for an ejection probability of p c j = 0.9 after six generations. 
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Fig. 3. — BH occupation fraction a s a function of the galactic veloc ity dispersion a v for the Monte Carlo merger model with v CBC = 2cr„ 
and the recoil velocity distribution of Schnittman & Buonanno (2007) (upper panel). The solid curves correspond to an initial occupation 
fraction fo = 0.99 and the dashed curves /o = 0.4, with /; the occupation fraction after an average of i mergers per galaxy. Both models 
converge to a single value of /oo, shown as a t hick black curve. In the lower p anel, we show the occupation fraction for an average kick 
velocity three times as large as that predicted by Schnittman & Buonanno (2007). In all cases, we set the merger fraction in each generation 
as fm = 0.66. 



